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CHAPTER  1: 
Introduction 


The  principle  of  acoustic  pressure  reciprocity,  originally  developed  by  Lord  Rayleigh  [1], 
is  a  simple  yet  powerful  statement  that  has  found  many  applications,  particularly  in  the  area 
of  propagation  modeling.  The  relationship  is  governed  by 

Pi(r2)  _  p2{ ri) 

p(r2)  p(ri) 

where  p(r  i)  is  the  density  at  point  r  i  and  p\(r2)  is  the  Green’s  function  at  point  r2  given  a 
source  at  point  r  For  the  acoustic  linear  wave  equation,  the  Green’s  function  is  related  to 
the  solution  of  the  acoustic  potential  function,  or  the  pressure  field,  due  to  omnidirectional 
point  sources  and  receivers. 

By  relating  the  pressure  received  at  one  point  in  an  acoustic  field  due  to  a  source  at  another 
point  to  the  pressure  received  when  source  and  receiver  positions  are  exchanged,  the  scalar 
reciprocity  relationship  allows  two-way  active  propagation  calculations  to  use  the  same 
path  for  the  return  echo  that  was  calculated  for  the  outgoing  waveform.  It  also  allows  one¬ 
way  passive  calculations  to  determine  the  pressure  received  from  multiple  source  locations 
by  computing  the  acoustic  field  for  a  source  at  the  receiver  location.  These  and  other 
applications  of  the  reciprocity  principle  greatly  reduce  the  computations  required  to  predict 
hydrophone  performance.  There  are  some  conditions  for  which  reciprocity  does  not  hold, 
i.e.,  in  the  presence  of  ocean  currents,  but  despite  those  exceptions  it  remains  a  useful  tool 
for  the  acoustician. 

While  hydrophones  are  only  capable  of  measuring  the  scalar  pressure  field,  vector  sensors 
are  able  to  measure  the  directional  components  of  acceleration,  velocity,  or  displacement. 
These  components  of  the  acoustic  field  contain  information  that  is  lost  in  the  scalar  pressure 
field,  and  exploiting  that  information  requires  the  ability  to  model  acoustic  vector  fields.  In 
order  to  predict  the  performance  of  a  vector  sensor,  one  must  be  able  to  calculate  the  vector 
field  components  received  by  that  sensor  for  any  source  position  in  an  arbitrary,  range- 
dependent  environment.  This  application  would  be  greatly  simplified  by  invoking  acoustic 
reciprocity,  but  the  equation  for  reciprocity  commonly  given  in  textbooks  holds  only  for 
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scalar  fields  such  as  pressure  and  velocity  potential.  The  equation  is  not  valid  for  velocity 
field  components,  which  are  proportional  to  the  gradient  of  the  potential  field.  Therefore, 
one  cannot  simply  calculate  the  vector  field  for  a  point  source  and  assume  those  vector 
quantities  are  what  would  be  received  by  a  vector  sensor  located  at  the  source  position. 

However,  as  Rayleigh  noted  in  regard  to  the  principle  of  reciprocity,  “apparent  exceptions 
depend  on  a  misunderstanding  of  the  principle  itself,”  and  the  failure  of  the  scalar  reci¬ 
procity  equation  to  correctly  give  the  directional  components  of  vector  fields  is  simply  a 
misapplication  of  reciprocity.  In  fact,  a  vector-scalar  reciprocity  equation  can  be  derived 
that  accurately  predicts  the  received  vector  fields  for  any  source  location  in  an  arbitrary, 
range-dependent  environment.  That  vector-scalar  reciprocity  relationship  can  be  incor¬ 
porated  into  propagation  models  and  used  just  as  the  scalar  reciprocity  relationship  has 
been. 

This  thesis  demonstrates  in  Chapter  2  that  the  scalar  reciprocity  relationship  works  for  pres¬ 
sure  fields  but  not  for  velocity  fields.  It  then  derives  the  vector-scalar  reciprocity  equation 
in  Chapter  3  and  shows  that  it  correctly  calculates  received  vector  fields  in  Chapter  4.  The 
conclusion  in  Chapter  5  summarizes  the  results  and  outlines  additional  ways  the  vector 
reciprocity  relationship  can  be  used. 
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CHAPTER  2: 

Limitations  of  the  Scalar  Reciprocity  Equation 


It  is  tempting  to  extend  the  reciprocal  method  for  calculating  received  pressure  fields  to 
models  that  calculate  directional  field  components  such  as  vertical  and  radial  displacement 
and  velocity.  However,  this  approach  breaks  down  because  while  the  reciprocity  relationship 
as  stated  does  hold  for  the  displacement  and  velocity  potentials ,  which  are  proportional  to 
acoustic  pressure,  it  does  not  hold  for  their  individual  components,  which  are  proportional  to 
the  pressure  gradient.  This  will  be  demonstrated  with  three  examples:  analytic  expressions 
for  normal  modes  in  a  range-independent  environment,  numerical  results  for  that  same 
environment,  and  numerical  results  for  a  range-dependent  environment.  The  text,  table,  and 
figures  in  Sections  2.1  and  2.2  are  adapted  from  [2] . 1 


2.1  Analytic  Expressions  for  Pressure  and  Particle 
Velocity  in  a  Range-Independent  Environment 

The  acoustic  field  of  a  point  source  in  a  range-independent  waveguide  has  a  closed-form 
solution  that  demonstrates  which  components  of  the  vector  field  obey  reciprocity.  The 
pressure  takes  the  form  [3] 


P(r,  =  ^m(Zs)^,„(z)H^\krmr),  (2.1) 

u  m=  1 

where  yYm(z)  is  the  normalized  mode  shape  function  for  the  mrh  normal  mode  propagating 
with  horizontal  wavenumber  knn.  This  function  depends  on  source  depth,  zs,  receiver 
depth,  z,  and  the  horizontal  distance  between  them,  r.  Exchanging  the  source  depth  with 

‘Reprinted,  with  permission,  from  T.  J.  Deal  and  K.  B.  Smith,  “Modeling  acoustic  vector  fields  for  inverse 
problems,”  in  Proc.  2016  IEEE/OES  China  Ocean  Acoustics ,  Jan.  2016.  This  publication  is  a  work  of  the 
U.S.  Government  as  defined  in  Title  17,  United  States  Code,  Section  101 .  Copyright  protection  is  not  available 
for  this  work  in  the  United  States.  IEEE  will  claim  and  protect  its  copyright  in  international  jurisdictions 
where  permission  from  IEEE  must  be  obtained  for  all  other  uses,  in  any  current  or  future  media,  including 
reprinting/republishing  this  material  for  advertising  or  promotional  purposes,  creating  new  collective  works, 
for  resale  or  redistribution  to  servers  or  lists,  or  reuse  of  any  copyrighted  component  of  this  work  in  other 
works. 
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the  receiver  depth  produces  the  same  result  for  p(r,zs),  satisfying  the  requirements  of 
reciprocity. 

Particle  velocity  can  be  calculated  from  pressure  by  computing  the  spatial  gradient  v  (r,  z)  = 
j^Vp(r,  z).  Its  radial  and  vertical  components  are  [4] 

J  oo 

i b(r,z)  =  - -  V  'Vm(Zs)'l'm(z)krmH<l1\krmr),  (2.2) 

1 ’z(r.z)  =  - — 1 — 7  V  'IV ( 3 ) ° ' l|'" i7)  tf(p r ) .  (23) 

4"Po  & 

Like  pressure,  the  radial  velocity  (;r  is  unchanged  when  source  and  receiver  depths  exchange, 
so  it  also  obeys  reciprocity  in  this  range -independent  environment.  However,  Equation  (2.3) 
has  different  dependence  on  the  source  and  receiver  depths.  Unless  the  mode  shape  function 
is  proportional  to  its  first  derivative,  vertical  velocity  does  not  obey  reciprocity. 

In  a  Pekeris  waveguide  with  a  fluid  layer  of  depth  h,  sound  speed  c,  and  density  p  bounded 
by  a  pressure  release  surface  and  semi-infinite  fluid  bottom  with  sound  speed  Cb  and  density 
Pb ,  the  mode  shape  function  is 


XI 'm(z)  =  An  sin  ( kzmz ) ,  o  <  z  <  h, 


(2.4) 


where  kzm  is  the  vertical  wavenumber  for  the  mth  mode  in  the  fluid  layer.  The  mode 
normalization  factor  Am  is  given  by  [5] 


2po 


,  sin  (2 kzmh)  .po  sin2  (kzmh) 

1 - - t -  +  J 


2  A  ->ii  h 


Pb  kbzmh 


(2.5) 


where  kbzm  is  the  vertical  wavenumber  for  the  mth  mode  in  the  bottom  halfspace.  The 
pressure  field  can  then  be  expressed  as 


1 p(r,  z)  =  ^  J]  A;n  sin  ( kzmzs )  sin  (kz 


iz)  HiQ1)(kr 


nr). 


(2.6) 


m= 1 
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The  mode  shape’s  vertical  derivative  is 


d%„(z) 

dz 


—  kzmAm  cos  ( kzmz )  • 


so  the  expressions  for  the  radial  and  vertical  velocity  become 

-1 


vr(r,z) 


4  ojpz 

r  0  m=  1 


^  A m  sin  (kzmZs)  sin  ( kzmz )  krmH ^  (krmv\ 


|  oo 

vz(r,  Z)  =  - - y  V  A“  sin  ( kzmzs )  cos  ( kzmz )  k^H^ (krmr) . 

ti 


(2.7) 


(2.8) 

(2.9) 


This  mode  shape  function  is  90  degrees  out  of  phase  with  its  derivative,  so  exchanging 
source  and  receiver  positions  will  produce  drastically  different  results  for  vertical  velocity, 
but  it  will  leave  the  radial  velocity  unchanged. 


2.2  Numerical  Solutions  to  the  Range-Independent 
Environment 

The  acoustic  vector  field  for  the  range-independent  environment  can  be  demonstrated  by 
implementing  the  analytic  equations  for  the  normal  modes  with  specific  values  for  the 
environment  and  source  parameters.  The  water  layer  is  200  m  deep  with  sound  speed  1500 
m/s  and  density  1.0  g/cm3.  The  bottom  has  sound  speed  1700  m/s  and  density  1.5  g/cm3. 
For  a  continuous  source  at  25  Hz,  there  will  be  seven  trapped  modes.  This  waveguide  also 
supports  higher  order  modes  and  a  branch  line  integral,  but  their  effect  is  limited  to  the 
region  near  the  source,  and  the  trapped  modes  dominate  at  ranges  greater  than  a  few  water 
depths,  so  they  are  omitted  here  for  clarity.  The  horizontal  and  vertical  wavenumbers  for 
the  trapped  modes  are  listed  in  Table  2.1. 

These  wavenumbers  are  used  with  Equations  (2.6),  (2.8),  and  (2.9)  to  generate  the  plots  for 
pressure,  radial  velocity,  and  vertical  velocity  for  both  a  fixed  source  and  a  fixed  receiver 
in  Figure  2.1.  For  the  fixed  source  cases,  the  source  is  at  a  depth  of  100  m  and  range 
of  0  km,  and  the  field  components  are  calculated  for  all  receiver  ranges  and  depths.  For 
the  fixed  receiver  cases,  the  receiver  is  at  a  depth  of  100  m  and  range  of  0  km,  and  the 
field  components  are  calculated  from  all  combinations  of  source  depths  and  ranges.  In 
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Table  2.1.  Horizontal  and  vertical  wavenumbers  of  trapped  modes  for  25  Hz 
source  in  example  Pekeris  waveguide.  Source:  [2], 


Mode  m 

krm 

k-z.m 

1 

0.1038  +  0.00001 

0.0137  +  O.OOOOi 

2 

0.1011  +  0.00001 

0.0275  +  O.OOOOi 

3 

0.0962  +  O.OOOOi 

0.0413  +  O.OOOOi 

4 

0.0893  +  0.00091 

0.0548  -  0.0015i 

5 

0.0774  +  0.0024i 

0.0706  -  0.0026i 

6 

0.0595  +  0.0045i 

0.0863  -0.003  li 

7 

0.0269  +  0.01271 

0.1021  -0.0033i 

this  way,  the  fixed  source  fields  represent  the  response  due  to  a  single  source  location,  as 
a  function  of  range  and  depth,  and  the  fixed  receiver  fields  are  a  function  of  source  range 
and  depth,  representing  the  response  from  a  mobile  source  location  as  observed  at  a  single, 
fixed  receiver. 

Computing  these  fields  using  the  analytic  expressions  is  straightforward  because  they  are 
formulated  as  a  function  of  r,  z,  and  zs,  and  any  two  of  these  parameters  may  be  varied 
to  create  a  two-dimensional  plot  with  the  third  parameter  held  constant.  However,  many 
propagation  models  only  compute  the  field  as  a  function  of  r  and  z  for  a  fixed  zs.  Using 
these  models  to  calculate  the  fixed  receiver  field  is  more  complicated.  A  source  is  placed 
at  a  point  (r,  zs)  in  the  waveguide.  The  geometry  is  then  shifted  horizontally  such  that  the 
source  is  located  at  a  range  of  zero  and  the  receiver  is  a  distance  r  away.  The  propagation 
model  is  run,  and  its  solutions  for  field  components  are  evaluated  at  the  receiver  location 
(r,  z),  where  receiver  depth  z  is  fixed.  The  pressure  and  velocity  components  measured  at 
this  fixed  receiver  position  are  plotted  at  the  source  location  (r,  zs).  This  process  is  repeated 
for  all  source  locations  until  the  entire  field  is  completed  with  the  required  range  and  depth 
resolution.  In  this  way,  the  fixed  receiver  plots  indicate  the  pressure  and  velocity  response 
due  to  a  source  located  at  each  point  in  the  field. 
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All  pressure  plots  are  in  decibels  of  transmission  loss  referenced  to  one  Pascal  at  one  meter 
for  a  monopole  with  source  strength  of  one  Pascal  at  one  meter.  This  provides  a  normalized 
transmission  loss  of  zero  dB  at  one  meter.  Velocity  plots  are  in  decibels  referenced  to  one 
meter  per  second  for  the  same  source  strength.  Contours  are  drawn  at  five  dB  increments. 


Range  (km)  Range  (km) 


(a)  Fixed  Source  Pressure 


(b)  Fixed  Receiver  Pressure 


(c)  Fixed  Source  Radial  Velocity 


(d)  Fixed  Receiver  Radial  Velocity 
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-180 

-200 

1-220 

-240 

Figure  2.1.  Vector  field  components  in  the  Pekeris  waveguide  calculated  from 
analytic  expressions  for  trapped  modes  using  scalar  reciprocity.  Pressure  is 
in  dB/ /IPa.  Velocity  is  in  dB//lm/s.  Transmission  loss  contours  are  drawn 
every  5  dB.  Adapted  from  [2], 


As  seen  in  Figure  2.1a  and  Figure  2.1b,  the  pressure  field  is  identical  when  source  and 
receiver  positions  are  exchanged.  The  same  is  true  for  the  radial  velocity  in  Figure  2.1c  and 
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Figure  2.  Id.  However,  the  vertical  velocity  plots  are  completely  different.  The  field  for  the 
fixed  source  has  vertical  velocity  nodes  around  40  and  120  meters  that  are  consistent  with 
range,  while  the  vertical  velocity  nodes  for  the  fixed  receiver  show  no  such  consistency  with 
range.  The  fields  at  zero  depth  are  also  significantly  different,  with  the  fixed  source  field 
having  maxima  for  receivers  at  the  surface  and  the  fixed  receiver  field  having  minima  for 
sources  at  the  surface. 

The  Monterey-Miami  Parabolic  Equation  (MMPE)  model  [6]  is  used  to  further  demonstrate 
the  differences  in  the  field  components  for  a  fixed  source  and  a  fixed  receiver.  It  correctly 
captures  the  effects  of  leaky  modes  and  the  lateral  wave,  and  it  also  supports  range-dependent 
calculations.  The  fixed  source  case  requires  a  single  run  of  MMPE,  while  the  fixed  receiver 
case  requires  one  run  for  each  source  range  and  depth  in  the  field.  Figure  2.2  shows  the 
pressure,  radial  velocity,  and  vertical  velocity  for  a  fixed  source  and  fixed  receiver  in  the 
example  Pekeris  waveguide. 

As  seen  in  Figure  2.2a  and  Figure  2.2b,  the  pressure  field  is  identical  when  source  and 
receiver  positions  are  exchanged,  even  in  the  near  field  where  the  lateral  wave  and  leaky 
modes  have  greater  influence.  Beyond  2.5  km  the  pressure  field  settles  into  the  far  field 
pattern  predicted  by  Equation  (2.6).  The  same  is  true  for  the  radial  velocity  in  Figure  2.2c 
and  Figure  2. 2d.  However,  the  vertical  velocity  plots  are  again  significantly  different.  This 
is  true  both  in  the  near  field  and  far  field.  In  all  cases,  the  results  from  MMPE  are  in 
excellent  agreement  with  the  analytic  solution  for  the  Pekeris  waveguide,  with  the  exception 
of  the  near-field  components  that  were  not  included  in  the  trapped  mode  solution. 

Figures  2. 3-2. 5  quantify  the  error  between  the  reciprocal  fields  calculated  with  MMPE  using 
a  monopole  source  at  the  receiver  location  and  the  reference  fields  calculated  by  moving 
the  monopole  source  throughout  the  water  column.  Magnitude  difference  is  measured  in 
decibels,  so  it  is  effectively  the  ratio  between  the  two  fields’  magnitudes  and  will  be  zero 
dB  when  they  are  equal.  Phase  error  is  the  phase  angle  of  the  reciprocal  field  minus  the 
phase  angle  of  the  reference  field  measured  in  degrees.  The  histograms  for  magnitude  and 
phase  error  use  1  dB  and  1  degree  bins  for  all  the  computed  range  and  depth  data  points  to 
display  the  distribution  of  errors.  The  histogram  bin  amplitudes  have  been  divided  by  the 
total  number  of  field  points  in  the  calculation  to  normalize  them.  If  the  error  was  exactly 
zero  everywhere,  the  histograms  would  have  a  single  bar  at  horizontal  coordinate  zero  and 
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amplitude  1.  Box  plots  quantify  the  statistics  of  these  distributions  by  displaying  the  mean, 
lower  quartile,  and  upper  quartile. 

In  addition  to  the  error  plots,  first  and  second  order  moments  are  calculated  over  all 
data  points  in  the  water  column.  For  phase  errors,  the  mean  and  standard  deviation  are 
computed  in  degrees.  For  magnitude  errors,  the  mean  and  standard  deviation  are  computed 
in  dimensionless  percent  error  instead  of  dB. 


(a)  Fixed  Source  Pressure 


(b)  Fixed  Receiver  Pressure 


(c)  Fixed  Source  Radial  Velocity 


n-140 

-160 

-180 

-200 
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(d)  Fixed  Receiver  Radial  Velocity 


Range  (km) 

(e)  Fixed  Source  Vertical  Velocity 
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(f)  Fixed  Receiver  Vertical  Velocity 
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Figure  2.2.  Vector  field  components  in  the  Pekeris  waveguide  calculated 
with  MMPE  using  scalar  reciprocity.  Pressure  is  in  dB//lPa.  Velocity  is 
in  dB//lm/s.  Transmission  loss  contours  are  drawn  every  5  dB.  Adapted 
from  [2] . 
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As  expected,  the  reciprocal  pressure  field  matches  well  with  the  reference  field.  The  error 
magnitude  mean  is  0.618%,  and  the  error  magnitude  standard  deviation  is  10.0%.  The 
phase  error  mean  is  -0.0501  degrees,  and  its  standard  deviation  is  3.28  degrees.  The  points 
in  the  field  with  the  greatest  magnitude  and  phase  errors  correspond  to  nulls  in  the  pressure 
field,  which  are  caused  by  destructive  interference.  At  these  points  the  pressure  magnitude 
is  40-60  dB  lower  than  points  nearby,  so  errors  in  the  field  have  little  impact  on  calculations 
such  as  transmission  loss. 


Range  (km)  Range  (km) 


(a)  Magnitude  Error 


(b)  Phase  Error 


Magnitude  Error  (dB) 


0.8 


Phase  Error  (deg) 


(c)  Magnitude  Error  Distribution 


(d)  Phase  Error  Distribution 


(e)  Magnitude  Error  Statistics 


Figure  2.3.  Computational  error  calculating  the  received  pressure  field  in 
the  Pekeris  waveguide  using  MMPE  and  scalar  reciprocity.  Magnitude  error 
contours  are  drawn  for  ±1,  ±2,  and  +3  dB. 
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For  the  range-independent  Pekeris  waveguide,  Equation  (2.8)  predicts  that  the  radial  velocity 
field  should  be  reciprocal,  and  that  is  supported  by  the  error  plots  in  Figure  2.4.  The 
magnitude  error  mean  is  0.321%,  and  the  standard  deviation  is  10.9%.  The  phase  error 
mean  is  0.0222  degrees,  and  the  standard  deviation  is  3.73  degrees.  These  numbers  are  in 
good  agreement  with  the  statistics  for  pressure  error. 


Range  (km)  Range  (km) 


(a)  Magnitude  Error 


(b)  Phase  Error 


Magnitude  Error  (dB) 


0.8 


Phase  Error  (deg) 


(c)  Magnitude  Error  Distribution 


(d)  Phase  Error  Distribution 
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(e)  Magnitude  Error  Statistics 
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(f)  Phase  Error  Statistics 


Figure  2.4.  Computational  error  calculating  the  received  radial  velocity  field 
in  the  Pekeris  waveguide  using  MMPE  and  scalar  reciprocity.  Magnitude 
error  contours  are  drawn  for  +1,  ±2,  and  ±3  dB. 
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The  reciprocal  calculation  for  vertical  velocity,  however,  is  a  poor  match  for  the  reference 
solution  at  almost  all  ranges  and  depths,  as  shown  in  Figure  2.5.  Note  that  the  vertical  axes 
on  the  histograms  in  Figures  2.5c  and  2.5d  and  box  plots  in  Figures  2.5e  and  2.5f  have  been 
expanded  compared  to  the  histograms  and  box  plots  for  pressure  and  radial  velocity  due  to 
the  large  error  spread  for  vertical  velocity.  The  magnitude  error  mean  is  105%,  and  the  high 
standard  deviation  of  945%  shows  that  the  errors  have  a  large  spread.  The  phase  error  mean 
is  low  at  0.657  degrees,  but  the  standard  deviation  is  99.6  degrees,  which  indicates  the  low 
mean  is  due  to  an  even  spread  of  large  phase  errors  rather  than  a  tight  group  of  low  phase 
errors.  The  one  exception  to  the  large  magnitude  and  phase  errors  lies  on  the  line  where  the 
source  depth  is  equal  to  the  receiver  depth.  It  is  only  at  this  depth  that  the  vertical  velocity 
given  by  Equation  (2.9)  is  the  same  when  source  and  receiver  positions  are  interchanged.  It 
will  be  shown  that  this  result  does  not  hold  in  range-dependent  environments. 

The  error  statistics  for  calculating  the  reciprocal  pressure  and  velocity  field  components 
are  summarized  in  Table  2.2.  Distilling  values  calculated  over  an  entire  field  into  a  single 
number  tends  to  obscure  the  details  of  how  well  the  calculations  do  or  do  not  match,  but 
in  this  instance,  the  standard  deviation  values  provide  good  insight  into  which  reciprocal 
calculations  perform  better. 


Table  2.2.  Vector  field  calculation  error  statistics  using  scalar  reciprocity  for 
a  monopole  source  in  the  Pekeris  waveguide 


Pressure 

Radial  Velocity 

Vertical  Velocity 

mean 

0.618 

0.321 

105 

Magnitude  (%) 

std 

10.0 

10.9 

945 

mean 

-0.0501 

0.0222 

0.657 

Phase  (deg) 

std 

3.28 

3.73 

99.6 
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Magnitude  Error  (dB)  Water  Column  Points  Depth  (m) 


(a)  Magnitude  Error 


(b)  Phase  Error 


0 


Magnitude  Error  (dB)  Phase  Error  (deg) 


(c)  Magnitude  Error  Distribution 


(d)  Phase  Error  Distribution 


(e)  Magnitude  Error  Statistics  (f)  Phase  Error  Statistics 

Figure  2.5.  Computational  error  calculating  the  received  vertical  velocity 
field  in  the  Pekeris  waveguide  using  MMPE  and  scalar  reciprocity.  Magnitude 
error  contours  are  drawn  for  +1,  ±2,  and  +3  dB. 
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2.3  Numerical  Solutions  to  the  Range-Dependent 
Environment 

The  shelf  break  test  case  from  the  SWAM  ’99  workshop  [7]  provides  varying  bathymetry 
and  sound  speed  that  is  useful  for  demonstrating  the  limitations  of  the  scalar  reciprocity 
equation  for  predicting  vector  fields  in  an  arbitrary  environment,  as  shown  in  Figure  2.6. 
The  water  has  constant  density  of  1.0  g/cm3,  and  the  bottom  is  homogeneous  with  sound 
speed  1700  m/s,  density  1.5  g/cm3,  compressional  attenuation  0.1  dB/d,  and  no  shear.  The 
range-dependent  bathymetry  has  depth 

f  100,  r  <  5000 

Zb  =  ,  (2.10) 

[^r,  5000  <  r  <  20000, 

where  Zb  and  r  are  in  meters.  For  depths  z  <  Zb,  the  water  sound  speed  in  m/s  is  defined  by 


c(r,z)  =  1450  +  4.6T(r,  z)  -  0.055T2(r,  z)  +  0.016z  (2.11) 


where  T(r,  z)  is  the  temperature  in  degrees  Celsius  defined  by 

j  -  /  „//inn_~\  \  \  2 


T(r,z)  =5+T0(r) 


sinh(^pl)' 


V 


265 


(2.12) 


with  the  surface  temperature  in  degrees  Celsius  defined  by 

7b(r)  =  15  -  5  (l  +  tanh(10  -  ^))  .  (2.13) 

This  SWAM’ 99  environment  was  intended  to  represent  a  canonical  shelf  break  environment. 
The  receiver  is  located  on  the  shelf  in  the  water  column  at  a  depth  of  35  m  and  range  of  0 
m.  The  goal  is  to  calculate  the  vector  field  received  by  that  sensor  for  a  25  Hz  monopole 
source  located  at  any  depth  in  the  water  column  for  ranges  from  0  km  to  20  km. 


The  vector  field  components  are  plotted  in  Figure  2.7.  On  the  right  are  the  reference  solutions 
calculated  by  placing  a  monopole  source  at  each  location  in  the  field  and  computing  the 
pressure  and  velocity  at  the  fixed  receiver  position  on  the  shelf.  On  the  left  are  the  pressure 
and  velocity  at  each  point  in  the  field  for  a  monopole  source  on  the  shelf. 
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Figure  2.6.  Range-dependent  geometry  showing  fixed  receiver  location, 
variable-range  source  depth,  bathymetry,  and  sound  speed  profile.  Sound 
speed  contours  are  drawn  every  5  m/s. 


Figure  2.8  shows  the  error  in  estimating  the  received  pressure  field  using  the  reciprocal 
method  with  a  monopole  source.  As  expected,  there  is  excellent  agreement  between  the 
reciprocal  method  and  the  reference  solution.  In  fact,  it  is  this  excellent  agreement  in  highly 
range-dependent  environments,  along  with  the  significant  computational  savings  achieved, 
which  makes  the  reciprocal  pressure  calculation  so  popular.  Pressure  error  magnitude  for 
all  data  points  in  the  water  column  has  a  mean  of  1.85%  and  standard  deviation  of  31.1%, 
and  the  pressure  phase  error  has  a  mean  of  -0.0228  degrees  and  standard  deviation  of  10.3 
degrees.  Further  inspection  of  the  error  surfaces  in  Figure  2.8  indicates  that  the  regions  of 
greatest  magnitude  and  phase  error  are  at  points  close  to  nulls  in  the  pressure  field  pattern. 
At  these  locations  of  destructive  interference,  the  pressure  magnitude  is  40-60  dB  less  than 
the  pressure  at  comparable  ranges,  so  errors  in  the  field  calculation  at  these  points  are  more 
easily  tolerated  in  most  applications. 

The  error  in  calculating  the  radial  velocity  using  the  reciprocal  method  with  a  monopole 
source  is  shown  in  Figure  2.9.  There  is  generally  good  agreement  between  the  reciprocal 
calculation  and  the  reference  solution  at  short  ranges  when  both  source  and  receiver  are 
on  the  shelf  and  the  environment  is  approximately  range-independent.  This  is  the  behavior 
predicted  by  Equation  (2.2)  and  seen  in  the  Pekeris  waveguide  example.  As  the  source 
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moves  beyond  the  range  where  the  bathymetry  and  sound  speed  profile  begin  to  vary,  the 
reciprocal  solution  begins  to  deviate  from  the  reference  solution.  For  all  water  column 
points,  the  magnitude  error  has  a  mean  of  6.18%  and  a  standard  deviation  of  51.9%,  and 
the  phase  error  has  a  mean  of  -0.405  degrees  and  a  standard  deviation  of  12.4  degrees.  As 
with  the  pressure  field  calculations,  the  locations  of  greatest  error  coincide  with  nulls  in  the 
radial  velocity  field  due  to  destructive  interference.  While  the  average  magnitude  and  phase 
errors  that  result  from  this  reciprocal  calculation  may  be  tolerable  for  many  applications, 


(a)  Fixed  Source  Pressure 


(b)  Fixed  Receiver  Pressure 


(d)  Fixed  Receiver  Radial  Velocity 


(c)  Fixed  Source  Radial  Velocity 


(e)  Fixed  Source  Vertical  Velocity  (f)  Fixed  Receiver  Vertical  Velocity 

Figure  2.7.  Vector  field  components  in  the  shelf  break  environment  calcu¬ 
lated  with  MMPE  using  scalar  reciprocity.  Pressure  is  in  dB//lPa.  Velocity 
is  in  dB//lm/s.  Transmission  loss  contours  are  drawn  every  5  dB. 
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the  trend  of  increasing  error  with  range  indicates  a  fundamental  flaw  in  using  a  monopole 
source  to  predict  received  radial  velocity. 

The  largest  errors  produced  by  this  misapplication  of  the  reciprocity  theorem  are  evident  in 
the  vertical  velocity  fields,  as  shown  in  Figure  2.10.  The  vertical  axes  on  the  histograms  in 
Figures  2.10c  and  2.10d  have  been  expanded  due  to  the  large  error  spread  and  consequent 
low  bin  amplitudes.  The  magnitude  and  phase  for  the  reciprocal  field  are  in  poor  agreement 
with  the  reference  solution  at  most  points  in  the  field.  This  wide  variation  leads  to  magnitude 
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(e)  Magnitude  Error  Statistics 


Figure  2.8.  Computational  error  calculating  the  received  pressure  field  in 
the  shelf  break  environment  using  MMPE  and  scalar  reciprocity.  Magnitude 
error  contours  are  drawn  for  +1,  ±2,  and  +3  dB. 
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and  phase  error  means  of  20.1%  and  0.261  degrees,  but  the  standard  deviations  are  much 
higher  at  1340%  and  98.7  degrees.  As  in  the  Pekeris  waveguide  example,  there  are  low 
magnitude  and  phase  errors  for  vertical  velocity  on  the  line  where  source  depth  equals 
receiver  depth  at  short  ranges  where  the  source  is  still  on  the  shelf  and  the  environment 
is  approximately  range-independent.  As  the  bathymetry  and  sound  speed  profile  begin  to 
change,  this  localized  region  of  low  error  vanishes,  so  it  is  not  a  feature  to  be  relied  upon  in 
range-dependent  environments.  Estimating  the  received  vertical  velocity  field  on  the  shelf 


Range  (km) 


(a)  Magnitude  Error 


(b)  Phase  Error 


20 

15 

10 
3 

JL  5 
H  ° 

|  -5 

CL 

-10 
-15 
-20 

(e)  Magnitude  Error  Statistics 

Figure  2.9.  Computational  error  calculating  the  received  radial  velocity  field 
in  the  shelf  break  environment  using  MMPE  and  scalar  reciprocity.  Magni¬ 
tude  error  contours  are  drawn  for  ±1,  ±2,  and  +3  dB. 


(f)  Phase  Error  Statistics 
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by  calculating  the  vertical  velocity  received  from  a  monopole  source  placed  on  the  shelf  is 
clearly  incorrect. 

The  summary  statistics  for  calculating  the  reciprocal  pressure  and  velocity  components  for 
the  range-dependent  environment  are  shown  in  Table  2.3.  The  accumulating  range  errors  in 
the  radial  velocity  calculation  are  hidden  in  the  single  numbers  for  error  mean  and  standard 
deviation,  but  the  vertical  velocity  calculation  errors  are  evident. 

These  examples  demonstrate  that  the  scalar  reciprocity  relationship,  as  stated  in  Equa¬ 
tion  1.1,  does  not  immediately  extend  to  cover  individual  components  of  the  vector  field. 
Instead,  a  new  reciprocity  relationship  must  be  derived  that  accounts  for  the  spatial  gradient 
of  the  pressure  field. 


Table  2.3.  Vector  field  calculation  error  statistics  using  scalar  reciprocity  for 
a  monopole  source  in  the  water  column  of  the  shelf  break  range-dependent 
environment 


Pressure 

Radial  Velocity 

Vertical  Velocity 

Magnitude  (%) 

mean 

1.85 

6.18 

20.1 

std 

31.1 

51.9 

1340 

Phase  (deg) 

mean 

-0.0579 

-0.334 

6.41 

std 

7.78 

11.4 

103 

19 


Magnitude  Error  (dB)  Water  Column  Points  Depth  (m) 


(e)  Magnitude  Error  Statistics  (f)  Phase  Error  Statistics 

Figure  2.10.  Computational  error  calculating  the  received  vertical  velocity 
field  in  the  shelf  break  environment  using  MMPE  and  scalar  reciprocity. 
Magnitude  error  contours  are  drawn  for  ±1,  ±2,  and  +3  dB. 
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CHAPTER  3: 

Vector-Scalar  Reciprocity  Equation  Derivation  and 

Applications 


This  derivation  follows  the  general  approach  of  Case  [8]  using  the  methods  in  [9].  All 
signals  have  an  implied  time  dependence  of  e~jajt  that  is  suppressed  for  brevity.  Bold  face 
indicates  vector  quantities. 


3.1  Reciprocity  for  General  Sources 

The  general  form  of  the  linearized  inhomogeneous  wave  equation  is  [10] 


P(  r) 


v  Vpi(r) 
p(  r) 


+  k2pi(r)  =  ja>Mi(r  -  n)  +  V  •  Fi(r  -  ri) 


(3.1) 


where  p\  (r)  is  the  pressure  at  r  given  a  source  at  ri  that  is  a  superposition  of  two  different 
types  of  sources.  The  first  source  is  an  acoustic  monopole  that  produces  mass  injection 
rate  per  unit  volume  M \ .  The  second  source  is  an  acoustic  dipole  that  produces  directional 
force  per  unit  volume  Fj.  Now  consider  two  solutions  to  Equation  (3.1):  p\  (r)  for  a  source 
at  ri  and  p2(r)  for  a  source  at  r2.  Multiply  the  first  equation  by  p2,  and  multiply  the 
second  equation  by  p\ ;  then  subtract  the  second  equation  from  the  first  and  integrate  over 
the  volume  of  interest, 


L 

-L 

-L 


Pi(  r)V 
P2(r) 


Vpt(r) 

P(  r) 


-Pi(r)V 


V/>2(r) 


dV 


P(  r) 
Pi(r) 
.  p(r) 


ja>M\{r  -  ri)  + 

P(  r) 


P(r) 

Fi(r-ri) 


jcoM2(r 


n)  ♦ 

p(  r) 


F2(r-r2) 


dV 


dV. 


(3.2) 


Applying  Green’s  Theorem  to  the  left  hand  side  of  Equation  (3.2)  converts  it  to  a  surface 
integral  that  vanishes  for  the  impedance  boundary  conditions  one  encounters  in  underwater 
acoustics  since  potentials  are  proportional  to  their  gradients.  Therefore,  the  integrals  on  the 
right  hand  side  must  be  equal. 
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This  is  the  form  of  the  general  reciprocity  relationship  for  pressure  and  velocity, 


XI 


po( r)  po( r) 

LAZjcoMl(r  -  n)  +  y-4^V  F!(r  -  n) 

P(r) 


P(r) 

L 


Pi(r) 

p(r) 


l'wM2(  r 


jy 


r2)  +  n«v 


P(r) 


F2(r  -  r2) 


jy.  (3.3) 


3.2  Reciprocity  for  Point  Sources 

Equation  (3.3)  can  be  further  simplified  for  point  sources.  Let  each  monopole  source  take  the 
form  M(r-ri)  =  M[  5(r-ri),  and  let  each  dipole  source  take  the  form  F(r-rj)  =  F i  c)  (  r — r  j ) 
where  «5(-)  is  the  three  dimensional  Dirac  delta  function.  Then  the  V  operator  can  be  changed 
from  F  to  p  via  integration  by  parts  and  the  Divergence  Theorem,  i.e., 


(V  •  A)  dV  =  - 


[  A  •  (Vp) 
Jv 

[  A  •  (yg) 
Jv 


dV  + 
dV  + 


dV 


When  g  =  />2(r)/p(r)  and  A  =  Fid(r  -  rj).  the  surface  integral  vanishes  since  the  delta 
function  is  zero  everywhere  on  the  surface,  and  the  volume  integral  becomes  the  dot  product 
of  the  force  and  the  gradient  of  the  pressure.  Equation  (3.3)  can  now  be  rewritten  as 


/[ 


jojM\5{v  -  rp 
P(  r) 


P2(r) 


XI 


Fi<5(r  -  r2)  •  V 

jcjMoSi r  -  r2) 
P(  r) 


(P2(r) 

\  P(r) 


dV  = 


Pt(r)  -  F2d(r  -  r2)  •  V 


Pi(r) 
P(  r) 


dV.  (3.4) 


Integration  over  the  entire  volume  extracts  only  those  values  where  the  delta  functions  are 
nonzero  at  the  source  locations  r2  and  r2,  resulting  in 


ja)Mip2(rj)  /P2(P)  \  =  jcoMm  (r2)  (piiXiL] 

p(ri)  1  \  p(ri)  /  p(r2)  2  \  p(r2)  / 


(3.5) 
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The  gradient  terms  in  this  equation  can  be  rewritten  as 


V 


=  V 


(  PC r)  \ 

\Po(l  +  J(r))/ 


(3.6) 


where  po  is  the  constant  equilibrium  density  and  s(r)  is  the  condensation.  In  [10]  the 
derivation  of  the  wave  equation  used  in  Equation  (3.1)  required  a  linearization  of  the 
equation  of  continuity  and  Euler’s  equation.  In  both  of  these  linearizations,  there  was  an 
assumption  that  s(r)  «;  1.  In  many  cases,  such  as  in  the  water  column  or  in  a  homogeneous 
layer,  that  assumption  is  still  valid,  so  the  gradient  term  can  be  reduced  to 


V 


Vp(  r) 

Po 


=  jojy(  r), 


(3.7) 


since  Vp(r)  =  ja>po\(r).  There  are  situations  where  this  assumption  does  not  hold,  such 
as  at  boundaries  between  layers  or  in  an  inhomogeneous  layer,  because  the  density  gradient 
is  not  exactly  zero.  In  these  cases  the  full  gradient  term  must  be  evaluated 


V 


Vp(  r) 
P(  r) 


+  V 


(3.8) 


Nevertheless,  returning  to  the  common  case  where  the  condensation  is  small  throughout 
the  water  column,  the  reciprocity  relationship  for  point  sources  can  be  stated  with  both  a 
pressure  component  and  a  velocity  component  as 

_  Fl  .  V2(ri)  =  M/-^  -  F2  •  vi(r2).  (3.9) 

P(ri)  p(Xi) 


3.3  Applications 

The  point  source  reciprocity  relationship  has  different  applications  depending  on  the  choice 
of  sources  used.  Two  monopole  sources  relate  pressure  to  pressure,  two  dipole  sources 
relate  velocity  to  velocity,  and  a  mix  of  monopole  and  dipole  sources  relates  pressure  to 
velocity. 
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3.3.1  Scalar  Reciprocity 

When  only  monopole  sources  of  equal  amplitude  are  considered,  then  M \  =  M2  and 
Fi  =  F2  =  0,  and  Equation  (3.9)  reduces  to 


Piirj)  _  p i(r2) 
p(ri)  p(r2)  ' 


(3.10) 


This  is  the  familiar  form  found  in  Equation  (1.1),  [1],  and  [9]  equating  the  pressure  at  point 
2  due  to  a  source  at  point  1  to  the  pressure  at  point  1  due  to  a  source  at  point  2. 


3.3.2  Vector  Reciprocity 

Similarly,  when  only  dipole  sources  are  considered,  M\  =  M2  =  0,  and  Equation  (3.9) 
becomes 

Fi  •  y2(n)  =  F2  •  vi(r2).  (3.11) 

This  is  the  reciprocal  relationship  between  velocity  and  force  components.  When  Equa¬ 
tion  (3.11)  is  expanded  in  Cartesian  coordinates  it  becomes 

F\xV2x<J\)  +  FiyV2y(r1)  +  F\zv2z(j\)  =  F2xv  i*(r2)  +  F2yv  i^(r2)  +  F2zv  iz(r2).  (3.12) 

This  equation  is  easier  to  use  when  the  forces  are  restricted  to  a  single  axis.  For  example, 
when  both  forces  act  only  in  the  x  direction,  Equation  (3.12)  reduces  to 

Fixv2x(n)  =  F2xv  i,(r2).  (3.13) 

When  the  amplitudes  of  the  force  components  are  equal,  F]x  =  F2x,  and  Equation  (3.13) 
states  that  the  x  velocity  received  at  point  2  due  to  a  force  in  the  x  direction  at  point  1  is 
equal  to  the  x  velocity  received  at  point  1  due  to  a  force  in  the  x  direction  at  point  2.  The 
same  applies  to  two  forces  in  the  y  direction  and  two  forces  in  the  z  direction. 

Equation  (3.12)  has  a  less  obvious  application  when  the  two  forces  act  in  orthogonal 
directions.  For  example,  when  the  first  force  acts  in  the  x  direction  and  the  second  force 
acts  in  the  y  direction,  Equation  (3.12)  becomes 

F\xV2X(j\)  =  F2yViy(r2).  (3.14) 
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When  the  force  components  are  equal,  F\x  =  F2y,  and  Equation  (3.14)  states  that  the  y 
velocity  received  at  point  2  due  to  a  force  in  the  x  direction  at  point  1  is  equal  to  the  x 
velocity  received  at  point  1  due  to  a  force  in  the  y  direction  at  point  2.  This  relationship 
holds  for  all  six  pairs  of  forces  in  orthogonal  directions. 

Equation  (3. 1 1)  is  useful  for  calculating  the  received  velocity  components  at  a  fixed  location 
due  to  a  dipole  source  at  any  range  and  depth  in  the  environment.  Suppose  the  dipole  source 
is  oriented  along  the  x  axis.  First  compute  the  x  velocity  everywhere  in  the  field  for  a  point 
force  with  unit  amplitude  at  the  desired  receiver  location  oriented  in  the  x  direction.  This 
quantity  is  represented  by  v(xX) ,  where  the  subscript  x  represents  the  field  component  that 
is  calculated,  and  the  superscript  (X)  represents  the  source  orientation  used  to  calculate  it. 
By  vector  reciprocity,  this  is  equal  to  the  x  velocity  received  at  the  desired  fixed  position 
for  a  force  in  the  x  direction  at  any  source  location  in  the  field.  Next  place  a  force  at  the 
desired  receiver  location  but  orient  it  in  the  y  direction.  Calculate  the  x  velocity  everywhere 
in  the  field,  vx  .  This  is  equal  to  the  y  velocity  received  at  the  desired  fixed  position  for 
a  force  in  the  x  direction  at  any  source  location  in  the  field.  Finally,  place  a  force  at  the 
receiver  location  oriented  along  the  z  axis.  Calculate  the  x  velocity  everywhere  in  the  field, 
vx  .  This  is  equal  to  the  z  velocity  received  at  the  desired  fixed  position  for  a  force  in  the  x 
direction  at  any  source  location  in  the  field.  Finally,  since  the  fields  were  calculated  with  a 
unit  amplitude  source,  multiply  each  quantity  by  the  dipole  source  amplitude  Fx. 

Since  any  force  can  be  decomposed  into  a  sum  of  components  in  each  orthogonal  direc¬ 
tion,  this  process  can  be  extended  to  arbitrary  source  orientations  via  superposition.  The 
reciprocal  velocity  field  must  be  calculated  for  each  orthogonal  source  direction.  This  can 
be  accomplished  with  three  propagation  model  runs  provided  all  three  velocity  components 
are  saved  for  each  model  run,  e.g.  vxX\  viX) ,  and  oiX>  for  the  source  in  the  x  direction. 
This  produces  three  x  velocities  ( vxx\  vx¥\  vxZ)),  three  y  velocities  (v(X) ,  v(¥\  v(Z)),  and 
three  z  velocities  (v(-X) .  v~¥\  v~Z)).  The  total  x  velocity  is  the  inner  product  of  the  three 
x  velocities  with  the  three  force  components,  the  total  y  velocity  is  the  inner  product  of 
the  three  y  velocities  with  the  three  force  components,  and  the  total  z  velocity  is  the  inner 
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product  of  the  three  z  velocities  with  the  three  force  components,  i.e., 


„  _  p  ,,(90  ,  17  .,09  ,  p  „(Z) 

Vx  —  r  yVx  +  t1  ZVX 

,,  _  p  „00  1  p  .,09  ,  p  ,.(Z) 

^1/  —  r XV  ij  +  r  yV  y  +  r  zVy 

vz  =  Fxv(zX)  +  Fyv(zY)  +  FzviZ).  (3.15) 


3.3.3  Vector-Scalar  Reciprocity 


The  final  reciprocal  relationship  arises  when  M 1  =  0  and  F2  =  0. 
tion  (3.9)  to 


-Ft  •  v2(ri)  =  M2 


Pi(r2) 
p(Xi)  ' 


This  reduces  Equa- 
(3.16) 


This  form  of  the  reciprocity  relationship  for  point  sources  permits  the  computational  short¬ 
cut  commonly  used  for  received  pressure  calculations  to  be  used  for  received  velocity 
calculations.  Normally,  a  monopole  source  is  placed  in  the  desired  receiver  location  and  the 
acoustic  field  is  calculated  everywhere  for  that  source.  This  pressure  is  then  scaled  by  the 
density  ratio  of  source  and  receiver  locations.  By  reciprocity,  this  quantity  is  the  pressure 
received  at  that  desired  fixed  position  for  a  source  located  at  each  point  in  the  field. 


Applying  this  method  to  Equation  (3.16)  is  a  two-step  process.  First,  a  dipole  source  is 
placed  in  a  particular  orientation  (aligned  with  the  x,  y,  or  z  axis)  at  the  desired  receiver 
location  and  the  pressure  field  is  calculated  everywhere  for  that  source.  The  pressure  at 
each  point  in  the  field  is  then  divided  by  the  density  and  negated  to  account  for  the  direction 
reversal.  Provided  the  sources  have  equivalent  magnitudes,  \Mo\  =  ||Fi||,  the  result  is  the 
amplitude  of  the  velocity  component  in  the  same  direction  that  the  dipole  was  oriented  in 
the  first  step.  To  obtain  the  other  received  velocity  components,  these  steps  are  repeated  for 
other  orientations  of  a  source  dipole  at  the  desired  receiver  location. 


3.3.4  Inverse  Vector  Fields 

When  combined  with  a  propagation  model  that  calculates  pressure  and  velocity  components 
of  the  vector  field  for  a  fixed  source,  Equations  (3 . 10),  (3 . 1 1 ),  and  (3.16)  enable  the  reciprocal 
computation  of  the  vector  field  for  a  fixed  receiver  from  the  solution  to  a  monopole  or  dipole 
source.  Let  the  model  calculate  the  received  pressure,  radial  velocity,  and  vertical  velocity 
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for  a  monopole  source  at  the  fixed  receiver  location,  p('M\  v[M  1 ,  and  vz.  Likewise,  the 
model  outputs  for  a  horizontal  dipole  source  are  pU!K  v(rH) ,  and  vzH\  and  the  model  outputs 
for  a  vertical  dipole  source  are  p(V\  v,<V\  and  viV>.  Using  these  model  outputs  properly 
scaled  for  units  of  Pascals  and  meters  per  second,  the  received  pressure  and  velocity  for  a 
fixed  receiver  and  monopole  source  are 


p  =  p(M> 


Vr  = 


1 


POC 

1 

POC 


,(ff) 


,09 


(3.17) 


the  received  pressure  and  velocity  for  a  fixed  receiver  and  horizontal  dipole  source  are 


P 

Vr 

VZ 


(POC)  V, 


( M ) 


AH) 


,09 


(3.18) 


and  the  received  pressure  and  velocity  for  a  fixed  receiver  and  vertical  dipole  source  are 


P  = 

Vr  = 


(p0c)  V 


(M) 


VZ 


(3.19) 


A  dipole  with  arbitrary  orientation  within  the  plane  of  calculation  can  be  parameterized 
by  its  rotation  from  vertical  by  an  angle  <p.  For  a  unit  magnitude  source,  the  horizontal 
and  vertical  components  of  this  dipole  are  equal  to  sin(^)  and  cos(</?),  respectively.  The 
received  pressure  and  velocity  for  a  fixed  receiver  and  arbitrarily-rotated  source  are 

p  =  (poc)  iy(M)  sin(^)  +  (p0c)  v\M)  cos(^) 
vr  =  vl-H^  sin(</?)  +  v[H)  cos (tp) 

vz  =  vj:V)  sin(</0  +  v(zV)  cos((f),  (3.20) 
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of  which  the  horizontal  and  vertical  dipole  are  special  cases  for  tp  =  |  and  <p  =  0.  These 
results  are  for  a  model  that  calculates  field  components  in  a  single  plane  with  one  vertical 
axis  z  and  one  horizontal  axis  r,  but  it  immediately  extends  to  a  model  that  calculates  a 
full  three-dimensional  field.  In  that  case  a  dipole  source  can  also  be  arbitrarily  rotated  by 
an  angle  ■&  in  the  horizontal  plane,  so  the  full  direction  cosines  needed  to  combine  model 
outputs  are  cos($)  sin(^),  sin($)  sin(^),  and  cos(cp). 
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CHAPTER  4: 

Demonstration  of  the  Vector-Scalar  Reciprocity 

Equation 


The  vector-scalar  reciprocity  relationship  in  Equation  (3.16)  can  be  demonstrated  analyt¬ 
ically  and  numerically  if  it  correctly  calculates  the  received  vector  fields  for  the  range- 
independent  and  range-dependent  environments  that  the  pressure  reciprocity  relationship 
failed  to  calculate  in  Chapter  2. 


4.1  Analytic  Expressions  for  Vector- Scalar  Reciprocity  in 
a  Range-Independent  Environment 

This  new  reciprocity  relationship  can  be  applied  to  the  range-independent  environment  from 
Section  2.1  by  calculating  analytic  expressions  for  the  pressure  field  created  by  horizontal 
and  vertical  dipoles.  Haug  et  al.  [11]  developed  an  expression  for  the  pressure  generated  by 
a  dipole  source  by  expanding  the  field  around  the  source  in  a  Taylor  series.  For  a  dipole  of 
source  strength  F\  -  ||Fi||  at  point  ri  rotated  </?  from  the  vertical  axis  and  f)  from  the  radial 
axis  in  the  horizontal  plane,  it  produces 

p(?2)  =  TT  'Yj  \-X¥m(z2)yi,m(zi)krmH[1) (knnr)  sin((p)  cos(P  -  6) 

J  m=0 

+  '¥m(Z2)'P'm(Zl)Hg\krmr)  cos(y>)]  .  (4.1) 

Here  r  =  \r2  -  r\ |  is  the  absolute  horizontal  distance  between  the  source  and  receiver,  which 
is  the  same  as  the  distance  r  in  Equations  (2. 1-2.3).  Equation  (4.1)  is  valid  for  any  mode 
shape  'T(z),  but  plugging  in  the  mode  shapes  for  the  Pekeris  waveguide  gives 


p(r2)  =  - y  ^  [~A;n  sin (kzmzi)  sin {k-mz\)knnH\l\knnr)  sin(^)  cos(tf  -  6) 


m= 0 


+ 


T  sin (kzmZ2)kzm  co^{kzmz\)H^\krmr)  cos(</?)]  .  (4.2) 
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When  the  dipole  is  oriented  horizontally,  sin(</?)  =  1  and  cos(</?)  =  0,  so  the  first  bracketed 
term  in  Equation  (4.2)  determines  the  pressure.  When  the  dipole  is  then  aligned  with  the 
plane  in  which  the  field  components  are  to  be  calculated,  cos  (7/  -6)  =  1,  and  the  pressure 
is 


P(  r2)  = 


-Fi 

4/ 


oo 

X  s  i  n  (7c  zm  Z2)sva{kzmz\)k 

rm  H?\k  rm ^  )  • 

m= 0 


(4.3) 


Similarly,  when  the  source  dipole  is  oriented  vertically,  cos(^)  =  1  and  sin(^)  =  0,  so  the 
pressure  field  is  determined  entirely  by  the  second  bracketed  term  in  Equation  (4.2), 


piX: 2)  =  —  'Yj  Am  Sin {kZmZ2)kzm  cos(kzmzi)H§\kr 

^  m= 0 


,r). 


(4.4) 


When  the  source  is  scaled  according  to  the  vector-scalar  reciprocity  relationship,  F\  = 
-M2/P0  =  -1/ po,  and  the  pressure  is  converted  to  velocity  through  division  by  jcopo,  these 
are  exactly  equal  to  the  velocity  expressions  in  Equations  (2.8)  and  (2.9).  Therefore,  the 
quantity  in  the  brackets  in  Equation  (4.2)  is  the  dot  product  of  the  dipole  source  with  the 
received  velocity  due  to  a  monopole  source,  as  predicted  in  Equation  (3.16). 


4.2  Numerical  Solutions  to  the  Range-Independent 
Environment 

Pressure  reciprocity  was  already  demonstrated  using  a  monopole  source  at  the  receiver 
location.  Now  the  reciprocal  field  for  radial  velocity  is  calculated  with  MMPE  using  a 
horizontal  dipole  source  at  the  receiver  location,  and  the  reciprocal  field  for  vertical  velocity 
is  calculated  with  MMPE  using  a  vertical  dipole  source.  The  method  used  to  modify  MMPE 
to  create  fields  for  dipole  sources  rather  than  its  standard  monopole  source  is  described  in 
the  appendix.  Since  the  received  field  to  be  calculated  here  is  for  a  monopole  source,  the 
reciprocal  velocity  fields  are  the  MMPE  pressure  output  for  each  of  the  directional  sources 
divided  by  the  characteristic  acoustic  impedance. 

Figure  4.1  shows  the  vector  field  components  measured  by  a  fixed  receiver  for  the  range- 
independent  Pekeris  waveguide  using  the  vector  reciprocity  relationship.  The  radial  velocity 
field  is  once  again  in  good  agreement  with  the  reference  solution.  This  is  borne  out  by 
the  error  analysis  in  Figure  4.2.  The  magnitude  error  for  radial  velocity  has  a  mean  of 
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1.31%  and  a  standard  deviation  of  9.67%.  The  phase  error  has  mean  -0.648  degrees  and 
standard  deviation  4.71  degrees.  These  numbers  are  comparable  to  the  statistics  for  the 
reciprocal  radial  velocity  calculated  with  a  monopole  source  in  Section  2.2,  but  that  is  what 
the  equations  predicted  for  a  range-independent  environment. 

As  compared  to  Figure  2.2,  the  reciprocal  vertical  velocity  field  in  the  left  column  is  a  much 
closer  match  to  the  reference  field  in  the  right  column.  Figure  4.3  shows  how  little  error 


(b)  Fixed  Receiver  Pressure 


(a)  Monopole  Source  Reciprocal  Pressure 
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(c)  Horizontal  Dipole  Source  Reciprocal  Radial  Velocity  (d)  Fixed  Receiver  Radial  Velocity 


(e)  Vertical  Dipole  Source  Reciprocal  Vertical  Velocity 


(f)  Fixed  Receiver  Vertical  Velocity 
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Figure  4.1.  Vector  field  components  in  the  Pekeris  waveguide  calculated  with 
MMPE  using  scalar  and  vector-scalar  reciprocity.  Pressure  is  in  dB//lPa. 
Velocity  is  in  dB//lm/s.  Transmission  loss  contours  are  drawn  every  5  dB. 
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there  is  between  solutions.  The  vertical  velocity  magnitude  error  mean  is  much  lower  at 
0.755%,  but  more  importantly,  the  standard  deviation  has  reduced  from  945%  to  14.0%. 
Similarly,  the  phase  error  mean  is  -0.398  degrees  and  the  standard  deviation  has  reduced 
from  99.6  degrees  to  5.42  degrees.  For  the  first  kilometer  of  range,  there  is  slightly  more 
variation  in  the  vertical  velocity  reciprocal  field  with  many  small  regions  of  fluctuation. 
Beyond  one  kilometer,  the  vertical  velocity  field  settles  into  its  far-held  pattern  and  the 
errors  are  restricted  to  locations  of  nulls  in  the  held  due  to  destructive  interference. 
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Figure  4.2.  Computational  error  calculating  the  received  radial  velocity  field 
in  the  Pekeris  waveguide  using  MMPE  and  vector-scalar  reciprocity.  Magni¬ 
tude  error  contours  are  drawn  for  +1,  ±2,  and  +3  dB. 
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Magnitude  Error  (dB)  Water  Column  Points  Depth  (m) 
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Figure  4.3.  Computational  error  calculating  the  received  vertical  velocity 
field  in  the  Pekeris  waveguide  using  MMPE  and  vector-scalar  reciprocity. 
Magnitude  error  contours  are  drawn  for  ±1,  ±2,  and  +3  dB. 
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Table  4. 1  lists  the  summary  statistics  for  error  in  calculating  the  received  vector  field  using 
the  vector  reciprocity  principle.  The  improvement  in  vertical  velocity  calculation  method 
is  evident  when  this  table  is  compared  to  Table  2.2. 


Table  4.1.  Vector  field  calculation  error  statistics  using  scalar  and  vector- 
scalar  reciprocity  for  a  monopole  source  in  a  Pekeris  waveguide 


Pressure 

Radial  Velocity 

Vertical  Velocity 

Magnitude  (%) 

mean 

0.618 

1.31 

0.733 

std 

10.0 

9.67 

14.0 

Phase  (deg) 

mean 

-0.0501 

-0.648 

-0.398 

std 

3.28 

4.71 

5.42 

These  results  come  at  a  substantial  computational  savings.  Generating  the  reference  field 
with  2-meter  depth  resolution  and  10-meter  range  resolution  required  40,100  model  runs  to 
calculate  the  received  pressure,  radial  velocity,  and  vertical  velocity  to  a  maximum  range 
of  4  km.  In  contrast,  the  same  resolution  was  achieved  in  3  model  runs  using  the  scalar  and 
vector-scalar  reciprocity  methods. 
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4.3  Numerical  Solutions  to  the  Range-Dependent 
Environment 

Returning  to  the  range -dependent  environment  from  Section  2.3,  the  baseline  reference 
fields  for  model  comparison  are  from  the  same  set  of  MMPE  runs  with  a  monopole  source 
at  depth  intervals  of  2  m  and  range  intervals  of  100  m.  The  reciprocal  pressure  is  calculated 
as  before  with  a  monopole  source,  but  the  reciprocal  vector  fields  are  now  calculated  using 
horizontal  and  vertical  dipole  sources.  The  resulting  fields  are  shown  in  Figure  4.4. 

Calculating  the  received  radial  velocity  field  using  the  reciprocity  relationship  for  a  horizon¬ 
tal  dipole  source  produces  an  excellent  match  to  the  reference  field,  as  shown  in  Figure  4.5. 
Unlike  the  incorrect  method  shown  in  Figure  2.9,  the  magnitude  and  phase  errors  do  not 
grow  in  range  but  remain  relatively  constant.  For  all  field  points  in  the  water  column,  the 
magnitude  error  mean  is  0.577%,  and  the  standard  deviation  is  29.2%.  The  phase  error 
mean  is  -0.187  degrees,  and  the  standard  deviation  is  10.6  degrees.  These  statistics  are 
comparable  to  the  statistics  for  the  reciprocal  pressure  calculation  error  in  Section  2.3. 

Even  more  pronounced  than  the  improvement  in  reciprocal  calculation  for  the  radial  velocity 
field  is  the  improvement  in  the  reciprocal  vertical  velocity  field  calculated  with  a  vertical 
dipole  source  at  the  receiver  location.  This  is  demonstrated  in  Figure  4.6.  In  stark  contrast 
to  the  incorrect  reciprocal  method  shown  in  Figure  2.10,  this  method  produces  a  close 
match  to  the  reference  solution  for  both  magnitude  and  phase.  The  magnitude  error  mean 
and  standard  deviation  for  all  field  points  in  the  water  column  are  5.44%  and  53.1%,  and 
the  phase  error  mean  and  standard  deviation  for  all  field  points  are  -0.505  degrees  and  10.0 
degrees.  As  with  the  radial  velocity,  these  error  statistics  are  comparable  to  the  reciprocal 
pressure  error  statistics. 
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(e)  Vertical  Dipole  Source  Reciprocal  Vertical  Velocity  (f)  Fixed  Receiver  Vertical  Velocity 

Figure  4.4.  Vector  field  components  in  the  shelf  break  environment  calcu¬ 
lated  with  MMPE  using  scalar  and  vector-scalar  reciprocity.  Pressure  is  in 
dB//lPa.  Velocity  is  in  dB//lm/s.  Transmission  loss  contours  are  drawn 
every  5  dB. 
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Figure  4.5.  Computational  error  calculating  the  received  radial  velocity  field 
in  the  shelf  break  environment  using  MMPE  and  vector-scalar  reciprocity. 
Magnitude  error  contours  are  drawn  for  ±1,  ±2,  and  +3  dB. 
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Figure  4.6.  Computational  error  calculating  the  received  vertical  velocity  field 
in  the  shelf  break  environment  using  MMPE  and  vector-scalar  reciprocity. 
Magnitude  error  contours  are  drawn  for  ±1,  ±2,  and  +3  dB. 
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The  statistics  for  the  reciprocal  vector  field  calculation  for  a  monopole  source  at  all  locations 
in  the  water  column  of  the  shelf  break  range-dependent  environment  are  summarized  in 
Table  4.2.  These  results  demonstrate  that  the  vector-scalar  reciprocity  relationship  is  valid  in 
arbitrary,  range-dependent  environments,  and  it  can  be  used  by  acoustic  propagation  models 
to  generate  accurate  vector  fields  at  a  fraction  of  the  computational  expense  required  for 
calculating  those  fields  point  by  point. 


Table  4.2.  Vector  field  calculation  error  statistics  using  scalar  and  vector- 
scalar  reciprocity  for  a  monopole  source  in  the  water  column  of  the  shelf 
break  range-dependent  environment 


Pressure 

Radial  Velocity 

Vertical  Velocity 

Magnitude  (%) 

mean 

1.85 

0.577 

5.44 

std 

31.1 

29.2 

53.1 

Phase  (deg) 

mean 

-0.0579 

-0.235 

0.543 

std 

7.78 

8.09 

3.03 

Once  again,  the  vector-scalar  reciprocity  method  provided  a  significant  computational  sav¬ 
ings.  The  reference  fields  for  the  shelf  break  environment  required  40,200  model  runs  to 
achieve  2-meter  depth  resolution  and  100-meter  range  resolution.  The  scalar  and  vector- 
scalar  reciprocity  methods  calculated  the  same  fields  in  3  model  runs  with  a  2-meter  depth 
resolution  and  10-meter  range  resolution.  Achieving  the  same  range  resolution  with  the 
moving  monopole  source  method  would  have  required  an  additional  360,000  model  runs. 
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CHAPTER  5: 
Conclusion 


The  fundamental  inability  of  the  classic  scalar  reciprocity  equation  to  accurately  predict 
vector  fields  stems  from  the  assumptions  in  its  derivation  of  monopole  sources  and  receivers. 
When  dipole  sources  are  included  in  the  reciprocity  derivation,  a  more  general  form  of  the 
equation  results.  This  equation  correctly  contains  the  scalar  reciprocity  relationship  between 
monopole  sources  and  pressure  fields,  the  vector  reciprocity  relationship  between  dipole 
sources  and  vector  fields,  and  the  vector-scalar  reciprocity  relationship  between  monopole 
sources,  dipole  sources,  scalar  fields,  and  vector  fields. 

The  scalar  reciprocity  equation  and  vector-scalar  reciprocity  equation  were  shown  to  accu¬ 
rately  predict  the  received  vector  field  created  by  a  monopole  source  in  a  range-independent 
Pekeris  waveguide  and  range-dependent  shelf  break  environment.  To  accurately  compute 
the  reciprocal  vector  fields,  the  parabolic  equation  model  MMPE  was  modified  to  produce 
vector  field  outputs  for  horizontal  and  vertical  dipole  sources. 

The  vector  reciprocity  equation  relating  vector  fields  to  dipole  sources  was  not  demonstrated, 
but  one  application  for  it  is  predicting  the  received  vector  field  created  by  wind-driven  surface 
noise.  This  noise  source  is  often  represented  by  a  large  number  of  vertical  dipoles  evenly 
distributed  on  a  horizontal  layer  less  than  a  wavelength  below  the  surface.  Propagating 
the  field  from  each  source  to  a  fixed  receiver  would  be  computationally  intensive  due  to 
the  large  number  of  sources  required,  but  by  invoking  the  vector  reciprocity  principle,  one 
propagation  model  run  with  a  dipole  source  at  the  receiver  location  could  predict  the  velocity 
received  from  the  surface  noise  sources  at  all  ranges. 

The  reciprocity  equation  derivation  can  also  be  repeated  with  higher  order  sources  such  as 
quadrupoles.  Although  like  dipoles,  quadrupoles  tend  to  be  inefficient  radiators  and  are 
not  commonly  used  in  underwater  acoustics,  their  inclusion  in  a  reciprocal  equation  could 
provide  some  mathematical  efficiency  for  calculating  other  quantities  of  interest  such  as 
field  circularity  or  higher-order  moments.  When  the  proper  relationship  between  source 
type  and  received  quantity  is  observed,  these  vector  and  vector-scalar  reciprocity  equations 
are  as  useful  as  the  scalar  reciprocity  equation. 
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APPENDIX:  Dipole  Starter  Fields  for  Parabolic 
Equation  Propagation  Models 


Acoustic  propagation  codes  that  use  the  parabolic  approximation  to  the  linearized  wave 
equation  are  commonly  implemented  as  one-way  marching  algorithms  that  require  an 
initial  condition,  or  starter  field.  One  way  to  define  a  starter  field  is  to  calculate  the  analytic 
solution  for  a  monopole  point  source  and  its  image  reflected  by  the  pressure  release  surface 
and  transform  that  to  the  wavenumber  domain  as  [6] 

/OO 

[6{z  -  zs)  -  6{z  +  zs)]e~jkzZdz  (A.l) 

OO 

=  -2 jaAM  sin (kzzs),  (A.2) 

where  a  is  a  normalization  factor,  Am  is  the  point  source  amplitude,  zs  is  the  source  depth, 
and  kz  is  the  vertical  wavenumber.  By  superposition,  the  starter  field  for  multiple  sources 
M  arranged  vertically  at  depths  zsi  is 

M- 1 

=  —2 joe  'Yj  Ai  sin (kzzsi).  (A.3) 

i= 0 

In  MMPE,  the  full  starter  field  includes  a  wavenumber-space  taper  function  (designed  to 
improve  phase  accuracy  at  higher  angles,  consistent  with  the  WKB  approximation)  and  a 
phase  function  that  accounts  for  grid  spacing, 

I Pm(kz)  =  -2 jyj ^  sin (kzZs)  |l  -  j  ejkz  *,  (A.4) 

where  the  normalization  term  and  source  amplitude  have  been  combined  into  the  term  with 
the  square  root.  To  take  advantage  of  the  vector  reciprocity  principle  in  MMPE,  additional 
starter  fields  must  be  defined  for  vertical  and  horizontal  dipoles.  Vecherin  et  al.  [12] 
prescribe  an  equivalent  source  method  that  decomposes  vertical  and  horizontal  dipole 
sources  into  vertical  monopole  arrays  with  element  amplitude  weights  such  that  the  array 
source  strength  is  equal  to  that  of  a  single  monopole.  Although  these  amplitude  weights 
and  element  spacings  can  be  implemented  with  MMPE’s  existing  vertical  array  source,  the 
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utility  of  dedicated  horizontal  and  vertical  dipole  sources  led  to  a  direct  implementation  of 
the  starter  field  in  the  vertical  wavenumber  domain  such  that  the  user  need  only  specify  the 
source  depth  and  frequency. 


A.l  Vertical  Dipole  Starter 

A  vertical  dipole  is  a  special  case  of  a  vertical  array  with  two  elements  180  degrees  out 
of  phase  with  each  other  and  separated  by  a  distance  2A z  much  smaller  than  a  wavelength 
( kAz  <  0.1).  The  desired  shape  of  the  acoustic  field  for  this  source  is  cos (cp)  for  angles  (p 
measured  from  the  vertical  axis.  The  acoustic  dipole  has  a  moment  equal  to 

D(cp)  =  -2jkAzAzcos(cp),  (A. 5) 


so  to  give  this  source  the  same  amplitude  as  a  monopole  with  constant  moment  D{(p)  =  Am, 
each  element  in  the  dipole  array  must  have  amplitude 


A- 


jAm 
2k  Az ' 


(A. 6) 


The  1 80  degree  phase  shift  is  accomplished  by  shifting  one  element  +90  degrees  and  shifting 
the  other  by  -90  degrees  with  an  additional  multiplication  of  ±j.  For  this  two-element  array, 
with  weights  +  jAz  at  depth  \zs  +  Az]  and  -jAz  at  depth  [zs  -  Az],  Equation  (A.4)  becomes 


<Pv(kz)  =  —2jyJ ^  [jAz  sin (kz[zs  +  Az])  -  jA-  sin (kz[zs  -  Az])]  |l  -  -A  j  ejk^i 

(A.7) 

When  this  equation  replaces  the  MMPE  monopole  starter  field,  the  model  outputs  pressure, 
radial  velocity,  and  vertical  velocity  for  a  vertical  dipole  source. 


A.2  Horizontal  Dipole  Starter 

The  horizontal  dipole  has  a  desired  shape  of  sin(0),  which  Vecherin  represents  as  a  sum  of 
cosines  for  angles  close  to  n/ 2  using  the  Taylor  series  approximation, 

sin(0)  *  1  -  -  cos2(0)  -  -  cos4(0).  (A.8) 
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Narrow-angle  parabolic  equation  codes  require  only  the  first  two  terms  in  Equation  (A. 8); 
wide-angle  parabolic  equation  codes  such  as  MMPE  require  all  three.  The  terms  in  Equa¬ 
tion  (A. 8)  can  be  represented  by  arrays  of  1,  3,  and  5  elements  vertically  spaced  by  2A z. 

Just  as  a  cos(0)  directivity  pattern  can  be  represented  by  a  superposition  of  two  monopoles 
180  degrees  out  of  phase  separated  by  a  small  distance,  cos2(0)  can  be  represented  by 
a  superposition  of  two  dipoles  180  degrees  out  of  phase  separated  by  a  small  distance. 
With  appropriate  spacing,  two  of  the  monopoles  can  occupy  the  same  position  and  their 
amplitudes  add,  which  means  only  three  monopoles  are  required  to  represent  cos2(0), 
where  the  amplitude  of  the  center  monopole  is  double  that  of  the  others  and  1 80  degrees 
out  of  phase.  This  process  extends  to  higher  order  terms,  and  in  general,  the  directivity 
pattern  cos"(0)  can  be  represented  with  a  vertical  array  of  n  -  1  monopoles.  The  repeated 
summation  of  equal  strength  monopoles  results  in  element  weights  that  follow  the  binomial 
coefficients  with  alternating  signs,  i.e., 

(-1  )kn' 

Ak  =  — — - — k  =  0  . . .  n,  (A.9) 

k\(n-k)\ 

where  k  is  the  zero-based  array  element  index.  As  with  the  vertical  dipole,  the  individual 
elements  in  this  array  must  also  be  scaled  by  the  amplitude  A"  =  to  produce  an 

overall  moment  amplitude  equal  to  the  monopole  moment  Am-  This  term  comes  directly 
from  the  Taylor  series  expansion.  Table  A.  1  lists  the  element  weights  required  to  produce 
each  cos "(0)  term  in  Equation  (A. 8)  including  both  Ak  and  A"  scaling  factors. 

Table  A.l.  Positions  and  weights  to  implement  cos n{(p)  directivity 


Depth  1 

COS  2(0) 

COS4(0) 

zs  -  4Az 

A4 

-  2A z 

A2 

-4A4 

zs  1 

-2A2 

6A4 

+  2A  z 

A2 

-4A4 

zs  +  4A  z 

A4 

To  approximate  the  sin(0)  directivity  pattern,  the  weights  for  each  element  in  the  vertical 
array  are  added  according  to  Equation  (A. 8).  This  results  in  a  five-element  array  with  the 


45 


weights  A0-4  shown  in  Table  A. 2.  As  with  the  vertical  dipole  starter,  this  approximation 
holds  for  ( kAz  <  0.1).  When  the  source  depth  is  close  to  a  boundary  in  terms  of  a 
wavelength,  this  representation  of  the  horizontal  dipole  as  an  extended  vertical  array  can 
create  problems  as  the  outermost  elements  can  cross  the  boundary.  While  adding  additional 
terms  to  the  Taylor  series  expansion  for  sin(0)  reduces  errors  at  higher  propagation  angles, 
it  comes  at  the  penalty  of  increasing  the  vertical  array  extent.  Therefore  the  best  choice  for 
array  size  is  five  elements  since  this  is  the  minimum  number  of  elements  required  to  meet 
the  wide-angle  requirements  of  MMPE. 

Table  A. 2.  Positions  and  weights  to  approximate  horizontal  dipole  source 

with  a  vertical  array  of  monopoles 


Depth 

~  sin(0) 

zs  ~  4A z 

A0  = 

_l/44 

zs  -  2A z 

Ai  = 

-|A2X  + 

I A4 

2  x 

zs 

Ao  = 

1  +  A2  - 

-  2  44 

4  x 

zs  +  2A  z 

a3  = 

-|A2  + 

1  /A 

2  x 

zs  +  4A  z 

a4  = 

-hAi 

For  implementation  in  MMPE,  each  of  these  array  elements  sums  with  its  image  source  to 
produce  the  starter  field 


$h(kz)  =  -2 j 


l go 
Ink 


[A0  sin (kz[zs  ~  4Az]) 


(A. 10) 


+A\  sin(fcj.[z,s  -  2Az])  +  A2  sin (kzzs)  +  A3  sin(A:z[z5  +  2Az]) 


A\4 


+A4  sin(fcz[z*  +  4Az])]  1  -  ejkz  f . 

k - 


With  this  starter  field,  MMPE  outputs  pressure,  radial  velocity,  and  vertical  velocity  for  a 
horizontal  dipole  source. 
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